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Abstract 

A search for the decays — >■ jj.^ fi~ fi'^ and B° — >■ is performed using data, corresponding to an integrated 

I luminosity of l.Ofb"^, collected with the LHCb detector in 2011. The number of candidates observed is consistent 

• • with the expected background and, assuming phase-space models of the decays, limits on the branching fractions are 

set: B{B° fi+^i-fi+fi-) < 1.6 (1.2) x 10"® and B{B° ti+^-n+fi-) < 6.6 (5.3) x 10"^ at 95% (90%) confidence level. 

In addition, limits are set in the context of a supersymmetric model which allows for the B^^-^ meson to decay into a 
> I scalar (S) and pseudoscalar particle (P), where S and P have masses of 2.5 GeV/c and 214.3 MeV/c, respectively, both 

resonances decay into fi^ fj,~ . The branching fraction limits for these decays are B{Bg — ^ SP) < 1.6 (1.2) x 10~* and 

Z3(B° -> SP) < 6.3 (5.1) X 10"'-* at 95% (90%) confidence level. 
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The decays B'^^-^ are strongly suppressed 

in the standard model (SM). Significant enhancements 
in the branching fractions can occur in beyond the SM 
theories [1] For example, in minimal supersymmet- 
ric models (MSSM) the decay is mediated via the new 
scalar 5", and pseudoscalar P sgoldstino particles, which 
both decay into iX^ n~ ■ The sgoldstinos can couple to 
SM particles via type-I couplings, where one sgoldstino 
couples to two SM fermions and via type-II couplings, 
where both S and P couple to two SM fermions in a 
four-prong vertex |T| . Such searches are also of particular 
interest since the HyperCP Collaboration found evidence 
of the decay E"*" — > pfi'^fi^, which is consistent with the 
decay E^ pP and P — > 1^-^ tJ-^ with a P particle mass of 
214.3 ±0.5 MeV/c^ ^Sj. The inclusion of charge conjugated 
processes is implied throughout this Letter. The Belle 
Collaboration has reported a search for a P particle in 
B° VP{-^ M+At") decays, where V is either a ii:*°(892) 
or a p°(770) meson [4]. The corresponding upper limits 
constrain the type-I couplings of sgoldstinos to the SM 
particles. The search presented here is sensitive to the 
decays B^^-^ — >■ SP in which both the S and P particle 

decay vertices are not significantly displaced from the B^^^ 
decay vertex and the P particle has a similar mass to that 
reported by the HyperCP collaboration. Such final states 
probe type-II couplings, which might dominate for general 
low energy supersymmetry breaking models and have not 
been probed by previous experiments. The present search 
is therefore complementary to the Belle study. Moreover, 
the four-muon final state is essentially background free and 
therefore ideally suited to search for new physics signatures. 
The present search is also sensitive to 5°^,^ — > 
decays which are not propagated through intermediate 
resonant structures. 

The dominant SM decay of a i? meson into a 
four-muon final state is B'^ — )■ J/-0(/)(lO2O), where both 
the J/ij} and the mesons decay into two muons 
[Fig. [l](a)]. In this Letter, this is referred to as 
the resonant decay mode. The branching fraction 
for — >■ J/i/) (— >■ (— > is calculated as 
the product of B{B^ ^ J/iPcj)), B{J/^ ^ n+ii-) and 
B{4>^ H+H') [5], resulting in a value of (2.3 ±0.9) x 10"^ 
The main SM nonresonant B^^-^ — ^ fi~ ^~ decay mode 
is i39s — > /i+/i^7 (— > where one muon pair is pro- 

duced via an electroweak penguin or box diagram and 
the other via a virtual photon [Fig. [l](b)]. The branching 
fraction of 5°^,^ — >■ (— >■ is expected to be less 

than 10^1" |6J. The diagram for the MSSM decay mode 
B° ->■ SP is shown in Fig.[l](c). 

This Letter presents a search for the decays 
^ /^^/^ M^A*^ using data corresponding to an in- 



tegrated luminosity of 1.0 fb of pp collision data at 
^ = 7 TeV collected with the LHCb detector in 2011. 
The resonant B^ J /ipcj) decay mode is removed in the 
signal selection and is used as a control channel to develop 
the selection criteria. The decay 5° — > J/ipK*'^, where 
J/ip ^ H'^ IjT and K*'^ — )■ iiT+Tr", is used as a normaliza- 
tion channel to measure the branching fractions of the 
^ A*^A* A^"*"/^" decays. 

The detector [7j is a single-arm forward spectrometer, 
covering the pseudorapidity range 2 < < 5, designed for 
the study of particles containing b or c quarks. The detec- 
tor includes a high precision tracking system consisting of a 
silicon-strip vertex detector surrounding the pp interaction 
region, a large-area silicon-strip detector located upstream 
of a dipole magnet, and three stations of silicon-strip detec- 
tors and straw drift tubes placed downstream. The magnet 
has a bending power of about 4 Tm. The combined track- 
ing system has a momentum resolution Ap/p that varies 
from 0.4% at momenta of 5 GeV/c to 0.6% at 100 GeV/c. 
The impact parameter (IP) resolution is 20 |J.m for tracks 
with high transverse momentum (pt)- Charged particles 
are identified using two ring-imaging Cerenkov detectors. 
Photon, electron and hadron candidates are identified by 
a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter. Muons are identified by a sys- 
tem composed of alternating layers of iron and multiwire 
proportional chambers. 

The trigger consists of a hardware stage, based on infor- 
mation from the calorimeter and muon systems, followed 
by a software stage that performs a full event reconstruc- 
tion. Events are selected by the single-muon, dimuon and 
generic 6-hadron triggers described in Ref. 8 . 

Simulated events are generated using Pythia 6.4 [S] 
configured with the parameters detailed in Ref. [lOj . 
Final state QED radiative corrections are included us- 
ing the Photos package [IT]. The EvtGen [15] and 
Geant4 [131 [13] packages are used to generate hadron de- 
cays and simulate interactions in the detector, respectively. 

Signal — > candidates are selected by ap- 

plying cuts on the final state muons and the reconstructed 
meson. Each final state muon candidate track is re- 
quired to be matched with hits in the muon system |15| . 
The muon candidates are required to have particle iden- 
tification (PID) criteria consistent with those of a muon 
and not those of a kaon or pion. This is determined by cal- 
culating an overall event likelihood for the distribution of 
Cerenkov photons detected by the ring-imaging Cerenkov 
system being consistent with a given particle hypothe- 
sis. To assess a particle hypotheses the difference between 
the logarithm of its likelihood and the pion hypothesis 
likelihood (DLL) is computed. Each muon candidate is re- 
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Figure 1: Feynman diagrams for the — !■ and — > ^i'^ fi^ fi'^ decays, via (a) the resonant — > J/i/^ip 

SM channel, (b) the nonresonant SM channel and (c) the supersymmetric channel. The latter is propagated by scalar S 
and pseudoscalar P sgoldstino sfermions. 



quired to have DLL{K - n)< and DLL(^ - 7r)> 0. The 
PID selection criteria yield a muon selection efficiency of 
78.5%; the corresponding efficiency for mis-identifying a 
pion (kaon) as a muon is 1.4% (< 0.1 %). All muon can- 
didates are required to have a track fit chi^ per degree 
of freedom of less than 5. Selection criteria are applied 
on the consistency of the muons to originate from a sec- 
ondary vertex rather than a primary vertex. The muons 
are each required to have the difference between the oi 
the primary vertex formed with and without the considered 
tracks, Xip, to be greater than 16. The i?^^,^ candidates are 
formed from two pairs of oppositely charged muons. The 
-6°^^ decay vertex is required to be less than 30 to en- 
sure that the four muons originate from a single vertex. In 
addition, the reconstructed B'^^^ meson is required to have 
a xip l6ss than 9 and hence be consistent with originating 
from a primary vertex. 

The B^^-j candidates are divided into two samples ac- 
cording to the invariant mass of the muon pairs. Sig- 
nal nonresonant candidates are required to have all four 
invariant mass combinations outside the respec- 
tive (j) and J/V' mass windows of 950 — 1090 MeV/ and 
3000 — 3200 MeV/c^. The four- muon mass resolution is 
estimated by the simulation to be around 19 MeV/c^ for 
both and B^ decay modes. The signal candidates 
are selected in a four-muon invariant mass window of 
±40 MeV/c^ around the world average 5°^^ mass 5 . Can- 
didate -Bg — )• J/ipcj) decays are used to optimize the selec- 
tion criteria described above. The B^ J /i'4> candidates 
are selected by requiring the invariant mass of one muon 
pair to be within the <f> mass window and that of the other 
pair to be within the J mass window. 

The selection criteria are chosen by applying initial cut 
values that select generic B meson decays. These values 
are then further optimized using B^ — > J/'ipcj) candidates. 



After applying the selection, seven B^ — >■ J/tpcf) candidates 
are observed in the signal sample. This is consistent with 
the expected 5° J/'^(j) yield, 5.5 ± 2.3, calculated by 
normalizing to the — > J/ipK*'^ decay mode. 

The dominant _B^^.j — > fi^ fi^ fi^ background is combi- 
natorial, where a candidate i?!^^^ vertex is constructed 
from four muons that did not originate from a single 
meson. Sources of peaking background are esti- 
mated to be negligible, the largest of these is due to 
B° ■iJj{2S) ^i+|J.-) K*° K+TT-) decays, which has 
an expected yield of 0.44 ± 0.06 events across the entire 
four-muon invariant mass range of 4776—5426 MeV/c^. 

To evaluate the combinatorial background, a single 
exponential probability distribution function (PDF) is 
used to fit the events in mass ranges of 4776—5220 and 
5426—5966 MeV/c^, where no signal is expected. Extrap- 
olating the PDF into the B'^ (Bg) signal window results 
in an expected background of O.SS^q^^ (0.30to;2o) events. 
Linear and double exponential fit models give consistent 
background yields. 

The branching fraction of the B^^^-^ — fi^ii^ fi^^^ de- 
cay is measured relative to that of the normaliza- 
tion channel i?° — )■ J/ipK*'^, this avoids uncertainties 
associated with the B^ production cross section and 
the integrated luminosity. This normalization channel 
has the same topology as the signal channel and two 
muons as final state particles. The S'-wave component 
from the nonresonant decay — ^ J /iI^K~^'k~ is removed 
from the present search. The branching fraction of 
Jjii) K*^ iiT+Tr-) is calculated as the 

product of B[B''^ — > J /ijjK*^), with the 5-wave component 
removed [B], B{K*^^ K+tt-) and ^ M+M") 0; 

the resulting branching fraction is (5.10 ± 0.52) x 10~^. 

The muon PID and kinematic selection criteria for the 
normalisation channel are identical to those applied to the 
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Figure 2: (color online). Invariant mass distribution of 
if"'"7r~/i+/^~ candidates. The i?" and B'^ signal distribu- 
tions are shown by short-dashed black and long-dashed 
red (gray) lines, respectively. The background shape is 
denoted by a long-dashed black line. The total fit result 
is shown as a solid blue (black) line. The inset shows the 
mass distribution centered around the i?? mass. 



signal channel. In addition to these criteria, the kaon is 
required to have its PID consistent with that of a kaon 
and not a pion, and vice versa for the pion. This re- 
moves i?" — >■ J/tpK*^ candidates where the kaon and pion 
mass hypotheses are exchanged. The K*'^ i.J/i') rneson 
is selected by applying a mass window of ±100 MeV/c^ 
(±50 MeV/c^) around the world average invariant mass [5]. 

There are four main sources of peaking backgrounds 
for — )■ J/ipK*^ decays. The first is the decay 
_B+ — > J/tjjK^ combined with a pion from elsewhere in 
the event, which is removed by applying a veto on candi- 
dates with a ^.^ fjT invariant mass within ±60 MeV/c^ 
of the world average B^ mass |S] • The second arises from 
B^ — >■ J /ip(f) decays, where the (/) meson decays to two 
kaons, one of which is misidentified as a pion. This back- 
ground is suppressed with a veto on candidates with a 
K~^iT~ invariant mass, where the pion is assigned a kaon 
mass hypothesis, within ±70 MeV/c^ of the world average 
(j) mass [S] . The third source of background is from the de- 
cay i?s — J/ipK*'^, which is included in the ■k~ fi^ n~ 
invariant mass fit described below. The last source of 
background is the 5'-wave component of the decay mode 
B^ — )• J/tjjK^TT~ , this is discussed later. 

The yield of the normalisation channel is determined by 
fitting the ^i^ invariant mass distribution with 

a combination of three PDFs, the B'^ J/^pK*'^ and the 
-Bg J/ipK*'^ signal PDFs consist of the sum of a Gaus- 
sian and a Crystal Ball function |17j . centred around the 



respective world average 5° and B^ masses and the back- 
ground PDF consists of a single exponential. The resulting 
fit is shown in Fig. [2} The B° J/ipK*° mass peak has 
a Gaussian standard deviation width of 15.9 ± 0.6 MeV/c^ 
and contains 31 800 ± 200 candidates. 

The branching fraction of the B'^^-^ — > decay 
is calculated using 



BiB] 
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where B{B'^ J/iIjK*°) is the branching fraction of 
the normalisation channel [5l[T6j and cb^^j /^k*'> and 
are the efficiencies for triggering, recon- 



structing and selecting the normalisation and signal chan- 
nel events, respectively. The efficiencies are calculated 
using simulated events and are cross-checked on data. 
The yields of the normalisation and signal channels are 
Nbo^j/^K'O and A^bo^j_^^+;,-^+^- , respectively. The rel- 
ative production fraction for B^ and B^_,^ mesons, fd(s)/ fd, 
is measured to be f^/fd = 0.256 ±0.020 for B^ decays [H] 
and taken as unity for i?" decays. The factor k accounts 
for the efficiency-corrected 5'-wave contribution to the nor- 
malisation channel yield; k is calculated to be 1.09 ± 0.09, 
using the technique described in Ref. (TH]. 

The PID components of the selection efficiencies are 
determined from data calibration samples of kaons, pions 
and muons. The kaon and pion samples are obtained from 
— > iir~7r+ decays, where the D° meson is produced 
via D'^tt^ decays. The muon sample is obtained 

from B~^ — J/ip{—i' fj,'^ fj,~)K~^ decays. The calibration 
samples are divided into bins of momentum, pseudorapidity, 
and the number of charged tracks in the event. This 
procedure corrects for differences between the kinematic 
and track multiplicity distributions of the simulated and 
the calibration event samples. 



Two models are used to simulate B' 







decays: (i) the phase space model, where the B^^-^ mass 
is fixed and the kinematics of the final state muons are 
distributed according to the available phase-space and 
(ii) the MSSM model, which describes the decay mode 
^(s) ~^ MSSM model the pseudoscalar particle 

P is a sgoldstino of mass 214.3 MeV/c^, consistent with 
results from the HyperCP experiment [3- The decay 
widths of S and P are set to 0.1 MeV/c^. The scalar 
sgoldstino S mass is set to 2.5 GeV/c^. If the mass of S 
is varied across the allowed phase space of the B'^^^ — >■ SP 
decay, the relative change in £3° ^^sp from the central 
values varies from —12.6% to +17.2%. 
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Table 1: Combined reconstruction and selection efficiencies 
of all the decay modes considered in the analysis. The 
uncertainties shown are statistical. 



Decay mode 


Model 


Efficiency [%] 




Phase space 


0.349 ±0.003 


B° SP 


MSSM 


0.361 ±0.003 


B° n+fi-fi+fi- 


Phase space 


0.359 ±0.003 


^ SP 


MSSM 


0.366 ± 0.003 


B" J/iiK*^ 


SM 


0.273 ±0.003 



Table 2: Systematic uncertainties on the branching frac- 
tions of i?!*^-) — >■ The combined systematic 
uncertainties are calculated by adding the individual com- 
ponents in quadrature. 



Source 


Systematic uncertainty [%] 


^(^0 ^ J/ijK*°) 


10.2 


S-w&ve correction 


8.3 


fd/ fs 


7.8 


Data-simulation differences 


5.2 


Trigger efficiency 


4.4 


PID selection efficiency 


4.1 


Simulation sample size 


1.3 


B° J/^K*° yield 


0.6 


Combined uncertainty 


17.2 


Combined uncertainty 


15.4 



The calculated efficiencies of all the simulated decay 
modes are shown in Table [U The total efficiencies of the 
MSSM models are comparable to those for the phase space 
models, indicating that the present search has approxi- 
mately the same sensitivity to new physics models, which 
feature low mass resonances, as to the phase space models. 

Systematic uncertainties enter into the calculation of 
the limits on B{B9.^ fi^ fj,~ fi~^ through the various 
elements of Eq. (fTl) . The largest uncertainty arises from 
the branching fraction of B'^ — )■ J/ipK*^, which is known 
to a precision of 10.2% [51 |TS]. An uncertainty arises 
due to the correction for the B° — >• J/iIjK~^tt~ S'-wave con- 
tribution. This is conservatively estimated to be 8.3%, 
which is the maximum relative change in k when it is 
calculated: (i) by using the angular acceptance from simu- 
lated events and (ii) by performing a fit with the physics 
parameters of the decay fixed and the angular acceptance 
parameterised, the coefficients of which are left free in 
the fit. An uncertainty of 7.8 % is introduced in the cal- 
culation of B{B^ — fj.^ ii~ fi'^ due to the uncertainty 
on fs/fd mi- The 4.4% systematic uncertainty associ- 
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Figure 3: Invariant mass distribution of nonresonant 
B{s) ~^ candidates. The soHd (dashed) black 

lines indicate the boundaries of the B^ {B^) signal window. 
The blue curve shows the model used to fit the mass side- 
bands and extract the expected number of combinatorial 
background events in the B^ and B^ signal regions. Only 
events in the region in which the line is solid have been 
considered in the fit. 



ated with the trigger efficiency is calculated as the relative 
difference between the data and simulation efficiencies of 
the trigger selection criteria applied to the normalisation 
channel. The efficiency in data is calculated using the 
method described in Rcf. 20J. Small differences are seen 
between the data and the simulated events for the track 
Xip distributions and the efficiency for reconstructing in- 
dividual tracks. The distributions of these quantities are 
corrected in the simulation to resemble the data using data- 
driven methods and the associated uncertainty is assessed 
by varying the magnitude and the configuration of the 
corrections. The relative systematic uncertainty assigned 
to the ratio of efficiencies, eB'>^j/t(>K*'> /^B°~>fj.+fj.- fi+ij.- ^ 
is calculated to be 5.2%. The 4.1% uncertainty on the 
PID selection efficiency is the maximum relative change 
in eBo^j/^K-'°/<^B'^^t,+ii-,j.+iJ.- that results from applying 
different binning schemes to the PID calibration samples. 
This uncertainty includes effects associated with muon, 
kaon and pion identification. The statistical uncertainty 
associated with the size of the simulated event samples 
is 1.3% for both B^ and B^ modes. The uncertainty on 
the B° J/ipK*" yield is 0.6%. Table [2I summarizes the 



systematic uncertainties for both the B{B^ — )• 
and B{B'^ — )■ 11^ ^~ fi~) branching fractions. The com- 
bined systematic uncertainty for the B^ and 5° modes is 
17.2% and 15.4%, respectively. The same uncertainties 
apply for ^(^0 ^ SP) and B{B° SP). 
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The nonresonant four-muoii invariant mass range and the 
background fit are shown in Fig. [3] One event is observed 
in the B'^ signal window and zero events are observed in 
the i?g window. These observations are consistent with 
the expected background yields. The CLs method pTll^ 
is used to set upper limits on the branching fractions. The 
95 % (90 %) confidence level limits for the nonresonant 
_B^^j decay modes in the phase space model 

are 

BiB^ fi+fi-fi+fi-) < 1.6 (1.2) X 10"^ 

B{B" ^1+ ^1- ^1+ ^i-) < 6.6 (5.3) X 10~^ 

The corresponding limits for the MSSM model with 
B°^^ SP and the mass of P{S) set to 214.3 MeV/c^ 
(2.5 GeV/c2), are 

B{B^ SP) < 1.6 (1.2) X 10-^ 

B{B° SP) < 6.3 (5.1) X 10"^ 

Varying the mass of S across the allowed phase 
space of the B°^^ ~> SP decay, from 211 MeV/c^ to 5065 
(5153) MeV/c^ for i?° (B^), results in a relative change in 
the 95% confidence level limit from —23% to +6% for 
both B° and B^ decay modes. 

In summary, a search for the decays i?"^,^ fi'^ fi^ jj.'^ j.i'' 
has been performed and first limits on the branching 
fractions for these decay modes have been set. These 
limits probe the upper regions of the parameter space 
of the B^^s^ SP decay, and, in particular set the first 
constraints on type-II couplings [1]. 
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